Lax Pairs and Darboux Transformations for 

Euler Equations 



Yanguang (Charles) Li * 
Department of Mathematics 

University of Missouri 
Columbia, MO 65211, USA 
e-mail: cli@math.missouri.edu 
phone: 573-884-0622 
fax: 573-882-1869 

Artyom V. Yurov 
Theoretical Physics Department 
Kaliningrad State University, Russia 
e-mail: yurov@freemail.ru 

February 1, 2008 



*This work is supported by the Guggenheim Fellowship. MSC: 76 34 35 37 



1 



Abstract 

In this article, we will report the recent developments on Lax pairs and Darboux 
transformations for Euler equations of inviscid fluids. 
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1 Introduction 

The governing equations for the incompressible viscous fluid flow are the Navier- Stokes 
equations. Turbulence occurs in the regime of high Reynolds number. By formally setting 
the Reynolds number equal to infinity, the Navier-Stokes equations reduce to the Euler 
equations of incompressible inviscid fluid flow. One may view the Navier-Stokes equations 
with large Reynolds number as a singular perturbation of the Euler equations. 

Results of T. Kato show that 2D Navier-Stokes equations are globally well-posed in 
C°([0, oo); H S (R 2 )), s > 2, and for any < T < oo, the mild solutions of the 2D Navier- 
Stokes equations approach those of the 2D Euler equations in C°([0, T]; H S (R 2 )) [10]. 3D 
Navier-Stokes equations are locally well-posed in C°([0, r]; H S (R 3 )), s > 5/2, and the mild 
solutions of the 3D Navier-Stokes equations approach those of the 3D Euler equations in 
C°([0, t]; H s (R 3 )), where r depends on the norms of the initial data and the external force 
[8] [9]. Extensive studies on the inviscid limit have been carried by J. Wu et al. [17] [6] [18] 
[4]. There is no doubt that mathematical study on Navier-Stokes (Euler) equations is one 
of the most important mathematical problems. In fact, Clay Mathematics Institute has 
posted the global well-posedness of 3D Navier-Stokes equations as one of the one million 
dollars problems. 

V. Arnold [1] realized that 2D Euler equations are a Hamiltonian system. Extensive 
studies on the symplectic structures of 2D Euler equations have been carried by J. Marsden, 
T. Ratiu et al. [14]. S. Friedlander and M. Vishik [7] [16] found a Lax pair for Euler 
equations written in the Lagrangian coordinates. 

Recently, the author [12] [13] found a Lax pair for 2D Euler equations written in the Eu- 
lerian coordinates. S. Childress [5] found a Lax pair for 3D Euler equations. In this article, 
a Darboux transformation for 2D Euler equations and their Lax pair is found. Darboux 
transformations are powerful tools in generating explicit representations for figure eight 
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structures in the phase spaces of Hamiltonian partial differential equations [11]. Of course, 
Darboux transformations were traditionally used for generating multi-soliton solutions to 
soliton equations [15]. 

Understanding the structures of solutions to Euler equations is of fundamental interest. 
Of particular interest is the question whether or not 3D Euler equations have finite time 
blow up solutions. T. Beale, T. Kato, and A. Majda derived a necessary condition [2]. 

Our hope is that the Lax pair and the Darboux transformation can be useful in in- 
vestigating finite time blow up solutions of 3D Euler equations and in establishing global 
well-posedness of 3D Navier-Stokes equations. 

The rest of this article is organized as follows: In section 2, we will present a formulation 
of Euler equations. In section 3, we will present a recent result on a Lax pair for 2D Euler 
equations. In section 4, we will present a recent result on a Lax pair for 3D Euler equations 
by Steve Childress. In section 5, we will present a recent result on a Darboux transformation 
for 2D Euler equations and their Lax pair. In section 6, a short conclusion is presented. 



2 Formulations of Euler Equtions 

The three-dimensional incompressible Euler equation can be written in vorticity form, 

d t n + (u-v)n-(n-v)u = o , (2.1) 

where u = (ui,u 2 ,u 3 ) is the velocity, Q = (Sli , Sl 2 , ^3) is the vorticity, V = (d x ,d y ,d z ), 
Q = V x u, and V • u — 0. u can be represented by Q for example through Biot-Savart law. 
When the fluid flow is of two dimensional, i.e. 

u = (u 1 (t,x,y),u 2 (t,x,y),0) , 

the vorticity Q is of the form Q = (0, 0, fls(t, x, y)). Dropping the subscript 3 of f^, the 2D 
Euler equation can be written in the form, 

d t n + {v,n} = o , (2.2) 

where the bracket { , } is defined as 

{f,9} = (dj)(d y g)-(d y f)(d x g) , 
where \& is the stream function given by, 

u = -d y ^ , v = d x ^ , 
and the relation between vorticity Q and stream function \1/ is, 

Q = d x v - d y u = . 
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(3-1) 



3 A Lax Pair for 2D Euler Equation 

Theorem 3.1 (Li, [12]) The Lax pair of the 2D Euler equation (2.2) is given as 

i Lip = \<p , 
\ dt<p + Atp = , 

where 

Lip = {{l,ip}, A(p = 
and X is a complex constant, and p is a complex-valued function. 

This theorem was announced in [12] without proof. Here we furnish a detailed proof. 

Proof: In order to solve the over-determined system (3.1), one needs a compatibility 
condition. To derive a compatibility condition, we first take a d t to the first equation, 

(d t L)<p + L{dt<p) = X(d t <p) = -XAp = -A(X(p) = —ALip , 

i.e. 

(d t L)(p — LA(p = —ALp , 

and finally 

(dtL) + [A,L^<p = 0, 
where [A, L] = AL — LA is the commutator. 

(d t L)p = {(d t n),p} , 

[A,L]ip = {%{n,<p}}-{n,{*M} 
{^,{n,p}} + {n,{p,^}} 



where by Jacobi identity, 



Finally 



0. Thus 

[A,L]<p = {{*,n}M- 

p = {d t n + {^,n},p} = o . 



(d t L) + [A,L\ 
Therefore, the 2D Euler equation 

d t n + {v,n} = o 

is a compatibility condition. □ 

Remark 3.1 A Darboux transformation for the above Lax pair will be presented in section 
5. Recently, Darboux transformations have been used in constructing explicit representa- 
tions of homoclinic structures [11]. 
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4 A Lax Pair for 3D Euler Equation 

Theorem 4.1 (Childress, [5]) The Lax pair of the 3D Euler equation (2.1) is given as 

{ Ly = \<p , 

\d t <p + A<p = 0, l4 ' ij 

where 

Lip = Q ■ Wtp — <p ■ Vf2 , Ap = u ■ Wtp — tp ■ Vw , 
A is a complex constant, and ip = (ipx, ip 2 , (ps) is a complex 3-vector valued function. 

This theorem was obtained by Steve Childress [5]. For a proof, see the proof of Proposition 
I. 

Remark 4.1 This Lax pair has the great potential that it can be very useful in investigating 
finite time blow up solutions of 3D Euler equation, and in establishing global well-posedness 
of 3D Navier- Stokes equations. 

Remark 4.2 When one does 2D reduction to the Lax pair (4-1), one gets L = 0. Therefore, 
the Lax pair (4-1) does not imply any Lax pair for 2D Euler equation. 

Remark 4.3 The most promising research following from this Lax pair should be along 
the direction of Darboux transformations, group symmetries etc.. In building the inverse 
scattering transform, it is crucial to have constant coefficient differential operators in (es- 
pecially the spatial part of) the Lax pair. 

Proposition 1 (Li, [13]) We consider the following Lax pair, 

i L<p = \<p, ( , 

\d tV + Ap = 0, 

where 

Lip = Q ■ Vy? — p ■ Vf2 + Di(p , Ap> = q ■ Vy? — <p ■ Vg + D 2 <p , 

X is a complex constant, p = (ipi, ip 2 , ipz) is a complex 3-vector valued function, q = 
(gi, q 2 , <?3) is a real 3-vector valued function, Dj = ■ V, (j = 1,2), and here = 
(ai\a 2 \a^) are real constant 3-vectors. The following equation instead of the 3D Euler 
equation, 

d t tt + (q ■ V)fi - (n ■ V)g + D 2 tt - D iq = , (4.3) 

is a compatibility condition of this Lax pair. A specialization of (4-3) is the following system 
of equations, 

d t n + (q- V)g = , 

D t q = D 2 tt . 
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A chance in building an inverse scattering transform [3] is that one can first build the 
inverse scattering transform for (4.2), and then take the limits 

a (»->0, (i = l,2) , q^u, 

to get results for 3D Euler equation. 

Proof of Proposition 1: Following the same argument as in the proof of Theorem 3.1, 
one has the compatibility condition, 



(d t L) + [A,L\ 



V? = . 



(dtL)<p = (d t Q) • Vy? - y? • V(d t Q) . 

[A,L]ip = q -V(fi -Vy? - (f -Vfi + Dup) 
—(Q • Vy? - <p ■ VVL + Diip) ■ Vg 
+D 2 (Q ■ Vy? - <p ■ Vfi + Dup) 
-tt ■ V(g • Vy? - ip • Vg + D 2 ip) 
+ {q ■ Vy? - ip • Vg + D 2 ip) • Vfi 
-D^q ■ Vcp - ip ■ Vq + D 2 ip) 
= l(q ■ V)fl] ■ Vy? + fl ■ [(g ■ V) Vy] (1) - [(g ■ V)y] ■ vn (2) 

-y. • [(q • V) Vfl] + (g-V)£>iy (q) - [(n-V)y]-Vg (4) 
+ [(y? • V)Q] • Vg - (Piy?)-Vg (?) + (L> 2 ft) • Vy? 
+ (^-V)D 2V 9 (6) - {D 2 ip)-Vn (7) - (y? • V)L> 2 ft 

+ ^ 2 ^iy (8) - [(n • v)g] • vy. - g-[(n-v)v y ] ( 1) 
+ [(n-v) y ]-vg (4) + <p ■ m ■ v) vg] - (n-v)fl 2 y (6) 

+ [(g-V)y]-Vfl (2) - • V)g] • Vfi + (D 2 y?)-V^ (7) 
— (D\q) ■ V<p - (g-V)% (3) + (%)-V 9 (5) 
+(cp-V)D 1 g- D 1 D 2 ip (B) . 
All the terms ^-), 1 < j < 8, cancel each other. Thus, one has 

[A, L]<p = j(g • V)fi — • V)g — -Dig + D 2 Q j ■ Vy? 

-if ■ V j(g • V)Q - (fi • V)g - Dig + D 2 nj . 

Therefore the following equation 

d t n + (g • V)fi - (J) • V)g + £> 2 fi - Dig = , 

is a compatibility condition. Since this equation involves two vaiables Q and g, one can 
choose the specialization, 

f d t ft + (g • V)fi - (fi • V)g = , 

\ Dig = d 2 q , 

to be a compatibility condition. □ 
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5 A Darboux Transformation for 2D Euler Equation 

Consider the Lax pair (3.1) at A = 0, i.e. 

{n, P } = o, (5.i) 

d t p+{q,p} = , (5.2) 
where we replaced the notation ip by p. 

Theorem 5.1 Let f = f(t,x,y) be any fixed solution to the system (5.1, 5.2), we define 
the Gauge transform Gj: 

P = G f p= ^-\p x - (d x In f)p\ , (5.3) 

and the transforms of the potentials fl and ^: 

y = ^ + F, tt=tt+AF, (5.4) 

where F is subject to the constraints 

{n,AF} = 0, {n + AF,F} = 0. (5.5) 

Then p solves the system (5.1, 5.2) at (Q, ty). Thus (5.3) and (5.4) form the Darboux 
transformation for the 2D Euler equation (2.2) and its Lax pair (5.1, 5.2). 

Remark 5.1 For KdV equation and many other soliton equations, the Gauge transform 
is of the form [15], 

P = Px~ {d x In f)p ■ 

In general, Gauge transform does not involve potentials. For 2D Euler equation, a potential 
factor ^- is needed. From (5.1), one has 

P^ = Py_ 

The Gauge transform (5.3) can be rewritten as 

_ = Px_ _ fx_P_ = Vy__ fy_P 
^x f ^y f 

The Lax pair (5.1, 5.2) has a symmetry, i.e. it is invariant under the transform (t,x,y) — > 
(—t,y,x). The form of the Gauge transform (5.3) resulted from the inclusion of the potential 
factor is consistent with this symmetry. 

Proposition 2 If F satisfies the constraints 

{n,AF} = 0, {n,F} = 0, (5.6) 

or the constraints 

{tt, AF} = , {AF, F} = , (5.7) 
then F satisfies the constraints (5.5). 
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Proof: Notice that {Q, AF} = implies that AF X = %^AF y . Thus, 
{Vt + AF, F} = (n x + ^AF y )F y -(n y + AF y )F x 

\ly 

= ^(Q y + AF y )F y -(Q y + AF y )F x 

y 



Similarly, one has 



Vly 



{Q + AF,F} = ^±^±{AF,F} 



Thus the claim in the proposition is true. □ 

Remark 5.2 For soliton equations, F can be represented in terms of f [15]. 
Proof of Theorem 5.1: One notices that using (5.1), (5.2) can be rewritten 



To prove the theorem, one needs to check the two equations, 

{(l,p}=0, 



■■■■ .'>•*> 



<-x 



First we check (5.9), 



leads to 



{n,p} = {n, pJ Q J fx } = o, 



where 



A 



nip nip 



(.Pxyf~\~Pxfy Pyfx pfxy)^xf (.Pxf Pfx){S^xyf~\~^xfy, 



Qjy 



B (Pxxf - Pfxx)^xf - (Pxf - Pfx){^xxf + fix fx) 

We will show that A — B. Using (5.1), one has 

Py\ (Vx 



Qy ) x V Q x / x 



which leads to 

Pxy xy Pxx Px^^x 

fijy fi« ^x ^x 
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Using (5.1) again, 
Similarly, 
Thus, 



Pxy Px^xy Pxx Px^i 



fxy fx^xy fxx fx^a 



Pxyf Pfxy (Pxf Pfx)^xy Pxxf Pfxx (Pxf Pfx)^x 



a 



That is, 



(Pxyf - Pfxy)ttxf _ {Pxf ~ Pfx)^xyf {j>xxf ~ Pfxx)^xf (Pxf - Pfx)®>xxf 



a. 



Uy *Ly 

Using (5.1) again, one has 

Pxfy Pyfx 



Using (5.12), (5.13), and (5.1), we have 

1 



fx Px r „ 

Px q (T" ~~ 



A = B 



f ^xPxx - f(ttxf)xPx + [fx(ttxf)x - fxx(Qxf)]p 



Thus equation (5.11) becomes 



If we let 

then (5.11) is satisfied. 
Next we check (5.10). 



A 



. 



nip 
= {n, AF} = o , 



pt = 



Pxf - Pfx 

~ <>,/ 



Using (5.8) for both p and /, one gets 

1 



Pt 



VIP 



\ 1 ij. / x \^l<x 



Pa 



~(Pxf -Pfx) 
{fPxx pfxx) 
HfPx-pfx) 



OJV ^ ^ ft2/2 



The right hand side of (5.10) is, 

Q x KJF J nj Q X + AF X 

PJx) nip Q X + AF X 
From (5.15) and (5.16), the coefficients of u fp xx — pf X x" leads to 

1 {fi,^} 1 {{I + AF, ^ + F} 

^^x f ^^x X f ^^x I P ^ 

and the coefficients of u fp x — pj x leads to 

1 ffrrw,^ -(^/U^ + AF,fr + F} 
OJV ^ A n x / 2 V l ' V, ft*/ 2 0, + AF, 

We will show that (5.17) implies (5.18). The left hand side of (5.18) is 

i s n x f{n,*} 
n x f\ n x ) x nip\ n x 



1 mM= ' R ftW t ) 



Thus (5.18) becomes 

i_ (Q n -(njun + AF,$ + F} 

fi 2 /2 i */;* ^ fi 2 /2 n x + AF x 

Therefore, (5.17) implies (5.18). From (5.17), we have 

(n x + AF x ){n, ^} = n x {n + af, * + f} , 

which leads to 

AF^ft, ^} = ft x [{ft, F} + {AF, F} + {AF, *}] . 

Notice that 

AF x {n, ^} - n x {AF, ^} = -v x {n, af} = o , 

by (5.14). Then (5.19) becomes 

Q x [{Q,F} + {AF,F}] = . 

Thus if 

{n + AF, F} = , 

(5.20) is fulfilled. □ 



10 



6 Conclusion 



We have briefly reported the most recent results on Lax pairs and a Darboux transformation 
for Euler equations of incompressible inviscid fluid flows. The most promising researches 
following from these Lax pairs should be along the directions of Darboux transformations, 
group symmetries, etc. These Lax pairs have the great potentials that they can be very use- 
ful in investigating finite time blow up solutions of 3D Euler equations, and in establishing 
global well-posedness of 3D Navier-Stokes equations. 
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